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ABSTRACT 



Context. Spectral analyses of hot, compact stars with NUIE (non-local thermodynamical equilibrium) model-atmosphere techniques 
allow the precise determination of photospheric parameters such as the effective temperature (T^ff), the surface gravity (log^), and the 
chemical composition. The derived photospheric metal abundances are crucial constraints for stellar evolutionary theory. 
Aims. Previous spectral analyses of the exciting star of the nebula A 35, BD-22°3467, were based on He-l-C-l-N-l-O+Si+Fe mod- 
els only. For our analysis, we use state-of-the-art fully metal-line blanketed NLTE model atmospheres that consider opacities of 23 
elements from hydrogen to nickel. We aim to identify all observed lines in the ultraviolet (UV) spectrum of BD-22°3467 and to 
determine the abundances of the respective species precisely. 

Methods. For the analysis of high-resolution and high-S/N (signal-to-noise) FUV (far ultraviolet, FUSE) and UV (HST/5r/5) obser- 
vations, we combined stellar-atmosphere models and interstellar line-absorption models to fully reproduce the entire observed UV 
spectrum. 

Results. The best agreement with the UV observation of BD-22°3467 is achieved at Teff = 80± 10 kK and log g = 7.2 ± 0.3. While 
T^ff of previous analyses is verified, logg is significantly lower We re-analyzed lines of silicon and iron (1/100 and about solar abun- 
dances, respectively) and for the first time in this star identified argon, chromium, manganese, cobalt, and nickel and determined 
abundances of 12, 70, 35, 150, and 5 times solar, respectively. Our results partially agree with predictions of diffusion models for DA- 
type white dwarfs. A combination of photospheric and interstellar line-absorption models reproduces more than 90% of the observed 
absorption features. The stellar mass is M x 0.48 Mq. 

Conclusions. BD-22°3467 may not have been massive enough to ascend the asymptotic giant branch and may have evolved directly 
from the extended horizontal branch to the white dwarf state. This would explain why it is not surrounded by a planetary nebula. 
However, the star, ionizes the ambient interstellar matter, mimicking a planetary nebula. 
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1. Introduction 

Abell 35 (g^n no ^ 1 2''53" '32':79, (Saooo = -22°52'22'.'55) was dis- 
covered by Abelll ( 1955L A 55 24) and classified as a planetary 



nebula (PNI. I Abelll (Il966h characterized it (A66 35, henceforth 
A 35, PNG303.6-I-40.0) as a homogeneous disk PN with an an- 
gular size between 636" and 938". Figure [1] shows an image of 
A 35 with very long exposure time. Its size is about 17' in east- 
west and 14' in north-south direction. The nebula shape appears 
to be quite atypical for a PN. The bow-shock structure is sur- 
rounded by a symmetric emission area that does not agree with 
the usual bi-polar or ellipsoidal morphologies. 

iGrewing & Bianchil (Il988h discovered a hot companion t o 



solved binary system ('Jacoby'l98l';'De Marco'2009'), composed 
of a W D star of spectral type DAO and a G8 III-IV - type com- 
panion dUnruh et al.ll200lt IStrassmeie3l2009l). The latter d omi- 
nates the spectrum at At 2800 A (iGrewing & Bianchilll988l) . 



the visible nucleus in the optical, (SAO 181201,'Jacobyl 



They suggested that it is an extremely hot DAO-typJ white 
dwarf (WD). A 35, BD-22°3467 (LWHya) is therefore a re- 



19811) . 



* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26666. 

** Based on observations made with the NASA-CNES-CSA Far 
Ultraviolet Spectroscopic Explorer. 
' spectral classification: characteristic H I and He ll absorption lines 



'Borkowsk i et al.l d 19901) explained the PN asymmetry 
(Fig.[T]) as the result of an interaction between the PN and the 
surrounding interstellar medium (ISM), the observed bow shock 
(the star is moving with v = 125 km/s) being located at the equi- 
librium sphere of ram pressure of the accelerated PN gas and the 
stellar wind ram pressure. 

A mean (HIPPARCOS, TYCHO-2, and UCAC-2 values) 
proper motion of /i„ = -60.9 5 ±1.5 mas/yr and ug - -14.63 + 
1.4mas/yr was measured bv iKerber et al.l (l2008h in agreern ent 
with HIPPARCOS values (HIP 62905. IVan LeeuwenI l2007h of 
fia = -60.91 ± 1.62mas/yr/itf = -13.41 + 1.24mas/yr. With 
its high proper motion, BD-22°3467 passes through the com- 
ple te visible nebula w ithin about 16000 years (Fig.[l])- Recently, 
Fr ew & Parked (1201 Oh claimed A 35 to be a "bow shock nebula 
in a photoionized Stromgren sphere in the ambient ISM ". Thi s 
scenario does not include a PN. Weidmann & Gamen (201 ih . 
however, classified A 35 as a binary PN (? + G8 IV) with a bc- 
CSPN (corresponding to binarity for the cool CSPN). 
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star. We begin with a brief description of the observed UV 
spectra and preparatory work (Sect. |2]l, followed by an intro- 
duction to our model atmospheres (Sect. O. Then we describe 
our spectral analysis (Sect. |4]i in detail and discuss the results in 
Sect.|5] 




Fig. 1. A 35: Narrow-band composite of Ho-, [Oiii], and [S ii] 
(in total 36 h exposure time ) by Dean Salman (Sh2-313, 
http://www.sharplesscatalog.com/). N is up, E is left, the FOV 
is 20' X 20'. The (red) arrow indicates the projected motion of 
BD-22°3467 during the last 10000 years. 



HIPPARCOS mea sured a parallax of 7.48 + 1.55 mas 



(iPerrvman et al.l[r997l D - 134;|^24Pc)- This was corrected later 



bv lGatti et alj (Il998h to D = 163''^^pc (this value was adopted 
bv lHerald & Bianchil2002l) . lGatti et al.l (Il998h determined a sep- 
aration of the hot and cool component of 13 - 28 AU at this 
distance. However, they adopted 160 pc as a minimum distance 
because previous estimates (see their Table 2) gave much higher 
value s (e.g. 360 + 80 pc from photometry measured by Jacobv 
ri98lh . fVan LeeuwenI (l2007i) presented a validation of the new 
HIPPARCOS Tsduction. The improved parallax (8.38+1.57 mas) 
is slightly larger with a similar error range. The corresponding 
distance of D = 1 19^^ l pc is now even smaller 

[Herald & Bianchil (2002) determined the atmospheric pa- 
rameters of both components and found Teff = 80 kK, log g-1 .1 
(cm/s^) for the WD, and reif = 5kK, log g = 3.5 for its com- 
panion star Their hot-component analysis of UV spectra 
{FUS^ HST/5774E and lU^ covering 905 A<^<3280A 
was performed with TLUSTY and SYNSPEC (iHubenvl Il988t 
lHubenv&Lan3ll992i 119951: iHubenv et al.lll994 . Thev deter- 
mined abundances for He, C, N, O, Si, and Fe (Tab.|4|. For the 
analysis of the cool companion they used LTE (local thermody- 
namical equilibrium) models (iKuruczill99ll) . 

State-of-the-art NLTE (Non-Local Thermodynamical 
Equilibrium) model atmospheres are fully metal-line blanketed 
and can consider all elements from hydrogen to nickel (RauchI 
[2003; Rauch et al. 2007). Therefore, we decided to re-analyze 
the available UV spectra of BD-22°3467 to determine abun- 
dances of hitherto neglected metals. The resulting abundance 
pattern can give clues to the evolutionary history of the DAO 



^ Far Ultraviolet Spectroscopic Explorer 

' Hubble Space Telescope / Space Telescope Imaging Spectrograph 



2. Observations and interstellar absorption 

Hot post-AGB (asymptotic giant branch) stars have their flux 
maximum in the UV wavelength range. Most of the exhibited 
metal lines are located there. Metal lines of successive ioniza- 
tion stages allow one to evaluate the ionization equilibrium of 
the respective species. This is a very sensitive indicator for T^ff. 
The error ranges of Tgff and abundance determinations can be 
reduced by the analysis of many spectral lines. 

The dominant ionization stages of the iron-group elements 
in the relevant T^ff and \ogg regime are v and vi (Fig.|2]i. Many 
strategic lines of these ions are located in the UV wavelength 
range. We therefore retrieved FUSE, HST/STIS, and lUE obser- 
vations from the MASTO archive (Tab. 1). 

Our spectral analysis is based on the high-resolution FUSE 
and STIS observat ions. The four STIS spec tra were co-added to 
improve the S/N. Herald & Bianchi (2002) analyzed the same 
FUSE spectrum but only one (O4GT02010) of the STIS spectra. 
The low-resolution lUI^ observations were used in addition to 
determine the interstellar reddening (^b-v, Fig- EJ as well as 
to verify the flux calibration. The STIS observation had to be 
scaled by a factor of 1 .995 to match the flux level of the lUE 
observation. This offset is possib ly due to bad centering of the 
exciting star during the exposure (iHerald & Bianchill2002h . 

2.1. Interstellar absorption and reddening 

The general procedure of determining the interstellar redden- 
ing is to normalize a theoretical spectrum as far as possible 
in the low-energy range because it is not significant there. In 
the case of BD-22°3467, the companion star dominates the ob- 
servation at wavelengths longer than ^ 2 800 A. Therefore, we 
constructed a theoretical, composite spectrum (with appropri- 
ate steflar radii, R'JR; = 7.43 ■ 10"^) using our final model 

spectrum (reff = 80kK, log^ = 7.2) and a Kurucz modeQ(reff = 
5000 K, logg = 3.5, solar abundances logZ = 0.0) for the hot 
and cool component, respectively. This spectrum was then nor- 
malized to the 2MASS K flux (Fig. |3]l. The best fit (matching 
all GALEX, HIPPARCOS, and 2MASS magnitudes and the UV 
spectra) was a chieved f or Es^ y = 0.02 ± 0.02, using the redden- 
ing law of Fitzpat ricfl (11999 *) with the stand ard j?v = 3.1. Our 
value agrees within error limits with those of iHerald & Bianchil 
(l2002i £b-v =0.04 ± 0.01), who used the Rayleigh- Jeans tail of 
the UV observations ( 1 000 A 2 ^ < 2 800 A) only. 

To determine the column density of interstellar neutral hy- 
drogen, He I, models with different nei were compared with 
Ly a. We found a best fit for hhi = 5.0 ± 1.5 x 10^° cm~^ ( Fig.|4l i, 
which agrees with the value found bv lHerald & Bianchil (l2002l) . 



2.2. Radial velocity of the hot component 

The STIS observation shows numerous sharp, isolated photo- 
spheric lines that are suitable to determine Vrad- On average (21 



International Ultraviolet Explorer 



' [http : //archive . st sci . edu/| 

* International Ultraviolet Explorer 

' [ftp : //ftp ■ stsci ■ edu/cdbs/grid/k93models/kp9S/l 
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Fig. 3. UV-optical-IR flux distribution of the BD-22°3467 bi- 
nary system, fitted by atmosphere models for both components. 
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Fig. 4. STIS observation around H i Ly a compared with our final 
model (Tab. |4]i at different «hi- The dashed line is the WD's 
atmospheric flux. 

lines of N V, O iv, O v, Ar v, Cr v, Cr vi, Fe v, Fe vi, and Ni v), we 
measured Vi-ad = - 14.4 ± 0.7 km sec"' . 

3. Model atmospheres and atomic data 

In this section, we briefly describe the atomic data and the pro- 
grams used for our analysis. Details about the calculation of the 
stellar atmosphere models can be found in Sect. l3.1l The usage 
of OWENS for the calculation of the ISM absorption is described 
in Sect.lT2l 

3.1. The photospheric model for the hot component 

Our model atmospheres (plane-parallel, chemically homoge- 
neous, in hydrostatic and radiative equilibrium) are calculated 
wkh TMAF^ the Tubingen NLTE model atmosphere package 
(fWerner etal] I2d03h . We calculated a small grid of models that 
include opacities of 23 elements from H - Ni (Tab.|2]i. 

|http : //astro . uni - tuebingen . de/- TMAP| 



H - Ar are represented by "classical" model atoms (iRauchI 

11997 1) taken from the Tubingen model atom database. 

For Ca - Ni, a statistical approach was applied to handle the 
large nu mber of atomic levels and fine transitions. We employed 
IrOnIc dRauch & Deetienll2003l: lRauch.2003,) to consider thou- 
sands of levels and millions of lines provi ded b y Kurucz (200^ 
and priv. comm.) and the Opacity Project (Seaton et al. 1994). 

Kurucz 's line lists are divided into so-called LIN (measured 
and theoretical lines) and POS (only measured, "good" wave- 
lengths) lists. For the model-atmosphere calculations, the LIN 
lists were used, to consider the total opacity properly. To cal- 
culate the emerging spectrum, POS lists were used to identify 
iron-group lines. Figure |5] demonstrates the difference between 
LIN and POS line lists. As an example for all species and their 
ions serves Crv (within the shown 1470 - 1500 A interval). At 
present, Kurucz provides 73 222 LIN and 249 POS lines of Cr v, 
459 LIN and 6 POS lines in this wavelength interval. We ad- 
justed the Cr abundance to [Cr] = -1-3.1 to reproduce the two 
strongest Crv POS lines, Crv /l/l 1482.76, 1489.71 A. At this 
abundance, most of the prominent Crv LIN lines appear dis- 
placed or simply too strong. An exception is Crv /1 1481.66 A, 
which matches the observation well. It can not be excluded that 
the Crv /i 1490.24 A LIN fines in the model is the observed 
1490.66 A line, for instance. However, reliable abundance de- 
terminations are only possible using unambiguously identified 
POS lines (Sect.l43b. 

In the framework of the Virtual Observatory (VCFl). all spec- 
tral energy distributions (SEDs, A - F^) calculated from our 
model-atmosphere grid are available in VO-compliant form from 
the registered VO service TheoSSA^ provided by the German 
Astmphysical Virtual Observatory (GAVCQ)- 

3.2. The ISM line-absorption model 

In the FUSE wavelength range (Fig. |6|l, numerous interstellar 
H2 and H i absorption lines hamper the analysis of the photo- 
spheric spectrum. We applied the program OWENS to calculate 
ISM line-absorption models. OWENS allows one to model dif- 
ferent clouds with distinct radial and turbulent velocities, tem- 
peratures, column densities, and chemical compositions. It fits 
Voigt profiles using a minimization. F or more information 
on OW ENS, see e.g. lH6brard et af] (l2002h or iHebrard & MoosI 
(l200l . 

Figure |6] shows a small section of the FUSE spectrum of 
BD-22°3467. Only the combination of the photospheric and 
ISM-model spectra reproduces most of its spectral features. Our 
ISM model includes lines of Hi, H2 (7 = - 9), HD, Cii- 
IV, Ni-iii, Nv, Oi, Ovi, Mgii, Aln, Sin-iv, Pii, Sii-iii, S vi. 
An -II, Mnii, Feu, and Niii. The radial velocity of the inter- 
stellar low-ionization (i - 11) atomic gas and H2 is found to be 
v'lad = 6 ± 2 km/ s. 

4. Analysis 

In this section, we begin with a verification of those atmospheric 
parameters that were already determined by Herald & Bianchil 
il2002h which may deviate because these authors used TLUSTY 
and considered He, C, N, O, Si, and Fe, while we employed 

' http : // astro . uni - tuebing en . de/- TMAD| 

http://www.ivoa.net 
" |http : //dc ■ g- VP ■ or g/theossa] 
fhttp : //www ■ g - VP ■ or g I 
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1480 1485 



Fig. 5. Section of the STIS spectrum (black line) compared with our TMAP model. The SED in the upper panel was calculated 
using Kurucz's POS data (the strongest lines are identified at top), the SED in the lower panel using LIN data. This section contains 
mainly Fe v - vi and Mn v - vi lines. The strongest Cr v LIN lines are identified at the bottom. 




930 935 940 945 950 955 960 965 970 975 980 985 990 995 



Fig. 6. Top: Section of the FUSE observation showing the necessity of a detailed modeling of the ISM. The dashed line is the 
final, pure photospheric model. The full (red) line shows a combined stellar and interstellar model. Blue marks at the bottom 
denote interstellar H2 absorption lines. Bottom: Comparison of (pure photospheric) TMAP models with log g -1.2 (thick, black) 
and log g-1 .1 (thin, red), normalized to the flux of the log g-7.2 model at 999 A for comparison. The full lines are calculated with 
the determined reddening of E'^^ - 0.04 (Sect.lTTT). the dashed lines (shifted in log/j by -0.4 for clarity) with Eg^^ = 0. 



TMAP with far more elements (Tab. |2|. This might result in a 
different at mospheric structure. We b egin with models with pa- 
rameters of iHerald & Bianchil (|2002|) and fine-tune them if nec- 
essary. 

4.1. Effective temperature and surface gravity 

IHerald & Bianchil (|2002|) determined reff = 80kK and 
log ^ = 7.7. TMAP test calculations showed that their T^ff 



value matches all ionization equilibria well, e.g. Niv/Nv, 
Oiv/Ov, and Mnv/Mnvi (Sect. l43b . A lower log^ = 7.2 
yields a significantly better fit agree of the theoretical 
He 11 /1 1640.42 A line profile with the observation in the 
line core, which is not matched at log g = 7.7 (Fig. |7]). For the 
line-profile calculations, we used the Stark line-broadening 
tables by Schoning & Butler ( 1989a b). The line wings are too 
broad at log g = 7.7 and somewhat too weak at log g = 7.2. The 
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same is true for Hen /1 1084.94 A (Fig. |7ll. The decrements 
of sp ectral line series such as the Hi Balmer (iRauch et alJ 
119981) . H I Lyman (transitions « - «' = 1 - 3 to 1 - 8), and He ii 
2-n' series (2-5 to 2-11) are strongly dependent on \ogg. 
The fit to the outer line wings of the latter two series in the 
FUSE spectrum is also very well reproduced at log g = 7.2. 
Figure |6] shows He ii 2-5 to 2-14 and the respective H i 1 - 3 
to 1-7 blends. A comparison of their theoretical profiles at 
logg = 7.2 and logg-1.7 with the observation shows that they 
are significantly too broad at log g = 7.7. This effect cannot 
be compensated for by a lower Zse-v (Fig. |6]l. The interstellar 
Hi absorption is not dominating the wings of these lines. The 
observed "shoulders" between the Hi Lyman lines and, thus, 
the Lyman decrement is well reproduced at log ^ = 7.2 and is 
definitely missed at log g -1.1. This is best visble between Lyf 
and LyS because the contamination of the FUSE observation 
with ISM lines is obviously weak, log § = 7.2 agrees with the 
value of Herald & Bianchi (2002) within realistic error limits of 
0.3 dex. Their assumed statistical error (logg = 1 .1^q\1) seems 
to be too optimistic. We finally adopted Teff = 80 ± iOkK and 
log g -1.2 + 0.3 for our analysis. 

4.2. Hydrogen and helium 

We adopted [He] = -0.3 ([X] den otes log [mass fraction / solar 
mass fraction] for species X) from lHerald & Bianchi (2002) and 
achieved a reasonable fit to the He ii fines (Fig.|7). 

4.3. Metal abundances 

In our TMAP model- atmosphere calculations, we considered 
opacities of all elements from H to Ni except for Li, Be, B, 
CI, and K. We convolved our synthetic spectra with Gaussians 
to simulate the respective instrument resolution. By comparison 
with the observed spectrum of BD-22°3467, we determined for 
the first time photospheric abundances of Ar, Cr, Mn, Co, and 
Ni. All identified lines of these species in the available obser- 
vations (Sect. |2]i were evaluated in our abundance analysis. No 
lines of F, Ne, Na, Mg, Al, Ca, Sc, Ti, and V were identified. 
For these, only upper limits could be derived by test models 
where the respective lines in the model emerge (at the abundance 
limit) from the noise in the observation. A detailed comparison 
of models with different abundances of individual metals showed 
that the typical abundance error is +0.3 dex. In the following, 
we briefly describe some measurements (ordered by increasing 
atomic weight). 

Carbon lines of stellar origin are not detectable in the spectrum 
of BD-22°3467. Narrow intersteflar C iv AA 1548.20, 1550.77 A 
lines are present (Fig. [8]) that virtually "bracket" the photo- 
spheric lines (vji^^j = -14.4km/s, Sect. l2.2l ). They have their 
origin in two different interstellar "clouds" with different ra- 
dial velocities and column densities (cloud 1: Vi-ad = -5.4km/s, 
log(«civ/(g/cm^)) = 16.8; cloud 2: v^d = -30.4km/s, 
log(nciv/(g/cm^)) = 14.6). The radial velocity of cloud 1 agrees 
with the velo city of th e low-ionization ISM gas (Sect. 13.21 ). 
Interestinglv. iSchneider et al.l ( 11983b determined Vjad - -6.6 + 
3.8 km/sec for the nebula. This even suggests a nebular origin of 
cloud 1 while cloud 2 with its significantly higher velocity lies 
on the line of sight toward us. 

It is worthwhile to note that no other interstellar two- 
component absorption-line feature is identified in the spectrum 
ofBD-22°3467. 
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Fig. 7. Spectra of our final model compared with the FUSE 
and STIS observations around He ii /1 1084.94 A (top, dashed: 
pure photospheric, full line: photospheric + ISM line absorp- 
tion) and He ii A 1640.42 A (bottom, dashed: log g - 1.1, full fine: 
log g = 7.2). The identified lines are marked. Most of the very 
strong absorptions in the top panel stem from H2. 



However, a photospheric component, even at low abun- 
dance, worsens the agreement between model and observa- 
tion in the shoulder between the interstellar lines (Fig. [8]l. 
Since the strongest Civ lines in the FUSE wavelength range 
(/1/1948, 1107, 1168A) become too strong for [C] > -2.9, we 
determine this as an upp er limit. This valu e is a bout four times 
lower than that found bv lHerald & Bianchil(l2002i) . 



Nitrogen Niv /1 1718.55 A and N v 1238.82, 1242.80A are 
prominent in the STIS observation. Our final model reproduces 
the strengths of these lines and the ionization equilibrium at 
[N] = -1.7 (Fig.© well. 

iHerald & Bianchil (|2002|) determined a much lower photo- 
spheric N abundance (< 10"^ times the solar value) because 
of strong line wings forming at higher abundances. The reason 
may be a different approximation by the quadratic Stark effect in 
TMAP and TLUSTY/SYNSPEC. 

For the line-absorption cross-sections TMAP calculates 



o"ii(v) = 



^/7Te^ /ij 
nieC Avd 



, H(a, v) 



(1) 
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Fig. 8. Civ AA 1548.20, 1550.77 A lines in the 5775 observation 
compared with our final model. Blue, thin: photospheric lines, 
red, thick: photospheric + interstellar lines, black, thin: inter- 
stellar lines. 




1243 



Fig. 9. Niv ^ 1718.55 A and Nv AA 1238.82, 1242.80 A lines in 
the STIS observation compared with our final model. 



TMAP uses fij of 7.8580 ■ 10"^ and 1.5716 ■ 10"', and classical 
damping constants T^d of 1.4401 ■ 10+'' and 1.4494 ■ 10+'' for 
N V 1242.80, 1238.82A, respectively. 

A comparison of calculated Fst values (Sahal-Brechot priv. 
comm.) has shown that rj"*''^^ / r|^^^a'-'^'"'='=''°' = 0.33 is a good 
approximation for the Nv resonance doublet within 10000 < 
r/K < 200000. We introduced this factor and adopted Tst = 
3.0 ■ rg™^'*^^ for our analysis. 



Oxygen is present in the spectra with the ionization stages iv 
and V. We can reproduce the ionization equilibrium and the 
line strengths. Figure [TOl shows Oiv AA 1338.61, 1343.51 A and 
O V 1 37 1 .29 A at [O] ^ -2.6. This agre e s well with the obser- 
vation and the value o f lHerald & Biaiichil dlool . 
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Fig.lO. Oiv 1338.61, 1343.51 A (left and middle) and 
O V /i 1371.29 A (right) lines in the 5775 observation compared 
with our final model. 



e and are electron charge and mass, c is the light velocity, Avp 
the Doppler width, /jj the oscillator strength for the transition i 
— » j, and 77 the Voigt function with 



Fluorine was not identified in the observation. We modeled 
Fv AA 1082.31, 1088.38 A and F vi ^ 1139.49 A. They would be 
visible in the observation at abundances higher than solar. 



a - 



and 



u = 



Av 
Avd' 



^^•^ Magnesium lines are weak and almost fade in the noise, 
even the strongest, Mgiv /1 1683.003 A. Figure |5] shows 
Mg IV A 1490.433 A, which is not detectable in the FUSE obser- 
vation at solar abundance. Thus, we can only determine an upper 

(3) abundance Hmit of [Mg] = 0.0. 



where Av = v - vp with the line fr equency vq. 

According to ICowlevI (Il970l [T97 1). the line-broadening due to 

the quadratic Stark effect is approximately given by 



- Stark 



= 5.5 X 10"^ 



z + 1 



(4) 



where is the electron density, T is the temperature, n"j^ is the 
effective principal quantum number of the upper level, and z is 
the effective charge seen by the active electron. Radiative (Fiad) 
and coUisional damping by electrons is considered so that 



r = F^H + F; 



Stark- 



Silicon exhibits only the Siiv AA 1393.76, 1402.77 A resonance 
doublet (Fig.fTTTi. other Si lines are too weak. A weak ISM con- 
tamination is present but we can determine [Si] = -2.1. This 
agrees well with the upper limit of lHerald & Bianchil(l2002l) . 



Argon is strongly enhanced ([Ar] = -t- 1 . 1 ) . Similar to the centra l 
star of the PN Sh 2-216, LS V -h46°21 (cf. lRauch etani2007h . 
Arvi AA 1283.96, 1303.89, 1307.38A are identified in the 5775 
observation (Fig. [T2l i and reproduced by our final model. Our 
models show a strong Arvii A 1063.55 A line but this is blended 
by interstellar H2 absorption and therefore not suitable for a Teff 
(5) determination via the ionization equilibrium. 
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strong change of the ionization balance toward Cr v is impossi- 
ble within the error limits, we have to consider these two lines 
as blends and not suited for a reliable abundance determination. 

Manganese exhibits several Mnv and vi lines in the 5775 
wavelength range, e.g. (the strongest) Mnv AA 1359.24, 
1382.88, 1432.84, 1440.31, 1443.31, 1452.88, 1457.46, 
1480.75, 1486.50 A and Mnvi AA 1272.44, 1285.10, 1333.87, 
1345.49, 1356.85, 1391.17, 1391.22, 1548.43 A (Fig. [HJ. All 
these lines are reproduced at [Mn] - +1.5. Since many lines 
of two successive ionization stages are suitable for spectral 
analysis, the Mn v - vi ionization equilibrium is one of the 
prime indicators for Teff for BD-22°3467. reff = 80kK is 
well-matched at log ^ = 7.2 (Sect. 14. It . 



Fig. 11. Si IV /l/l 1393.76, 1402.77 A in the 5775 observation 
compared with our final model. 

A problem are lines like Arv /t 1371.92 A (Fig. [T0| which 
are much stronger in our models than observed. Since many 
ionization equilibria of different species such as Mn v - vi 
(see below) are well-matched, a higher Teff that would shift 
the ionization balance from Arv to Arvi can be excluded. 
We checked the possibility that the oscillator strengths of 
Arv lines that we use in TMAP (from the Opacity Project, 
ISeaton et all 1 19941) are too strong. The Arv triplet 3p^ ^S" 
- 3p4 3p (Arv ^^1341.57, 1350.39, 1371.92A) is found in 
Opacity Project data (far off) at 1455.77 A (without fine-structure 
splitting) with / - 4.2333 ■ 10"^. From this value, we calculated 
/-values of 4.7037 ■ lO^^, 1.4111 ■ lO^^, and 2.3518 ■ 10"^ for 
the thre e fine-structure co mponents. A comparison with recent 
worko f Taval et al.l(l2009l f,.7, oo = 2.257-10-2 and 2.834-10-2) 
shows good agreement with the Opacity Project. Thus, the reason 
for the disagreement of the Ar v lines remains unknown. 




1283 1284 



1304 1305 

x/A 



1306 1307 1308 



Fig. 12. Section of the STIS spectrum around 
Arvii^l283.96, 1303.89, 1307.38 A compared with our 
final, composite photospheric and ISM spectrum. 



Chromium is about 70 times enhanced ([Cr] - +1.8). In 
Fig. [Uli, Cr VI AA 1417.66, 1455.28 A are shown that are well- 
reproduced by our final model. Crv 1482.76, 1489.71 A 
(Fig. |5]) would require a much higher [Cr] - +3.1 (Sect. 13.1b . 
Since T^ff and log g are well-determined (Sect. 14. Il l and, thus, a 




1424 



1425 



1426 



1427 



1428 



0.7 



0.6 



0.5 



0.4 



>>>>> >>>>>>>>>> >>> 

— O — Q) — OcD — O OJ — O — O CD ^ — ^ 




Fig. 14. Sections of the STIS observation containing (in addition 
to others) lines of Cr vi (top), Co (middle), and Ni (bottom) com- 
pared with our final model. 



7 



M. Ziegler et al.: BD-22°3467, a DAO-type star exciting the nebula Abell 35 




Iron h as many lines of Fe v and vi in the UV spectrum. In con- 
trast to iHerald & Bianchil (f2002:). who determined a 0.1 - 0.5 
times solar Fe abundance, we find it approximately solar (e.g. 
Fig. |5]l. There are, however, some Fe lines that are stronger in 
our model than observed (see some Fe vi lines in Fig. [14] bottom 
pa nel). Within the error ran ge, we find the same Fe abundance 
as lHerald & Bianchil (l2002h . 

Cobalt is strongly enhanced. Cow AA 1331.50, 1349.60A indi- 
cate [Co] = -H2.2 (Fig.IH). 

Nickel is also overabundant. Numerous Niv and vi lines are 
found in the UV spectrum. They are all reproduced at [Ni] - 
+0.7 (e.g. Figs.|5][Il. 

4.4. Mass, luminosity, radius, and distance 

We determined the stellar mass and luminosity from a com- 
parison with evolutionary tracks for hydrogen-rich post-AGB 
stars (Fig.[T5j Miller Bertolami, priv. com m.) and for po s t-EHB 
(extended horizontal branch) stars from iDorman et alJ (1 19931 
Fig- [13 ■ The results are summarized in Tab. [3] 

The stellar radius follows from R - ^jGM I g, where G 
is the gravitational constant. We determined the spectroscopic 
distance D - 361^[3^pc (Table|3]l following IWassermann et al] 
(12010) . It is derived from 

D/pc = 6.6175 ■ 10"'' y/Fy/F„b.- M/Ms ■ lO-^°ss 

with the astrophysical flux Fy at the surface of the star as 
predicted by the model spectrum, and Fobs is the observed, 
extinction-corrected flux (we used ^^1513^ to match the center 
of the GALEX FUV band). D is about three times larger than the 
improved ///PPA/;C05 distance (D = 119!j9pc, Sect. [B but it 
agrees with the photometric distance (360 ± 80 pc. Ilacobvl 1 98 iT) 
for the cool component. 




5.3 5.2 5.1 5.0 4.9 4.8 4.7 

log 7-,,, / K 



Fig. 15. Location of BD-22°3467 in the log Teff-logg di- 
agram c ompared with other DAO-type WDs taken from 
Gianni nas et alj ( l20IOl all symbols). Three of these (green 
squares) are used to compare their metal abundances with those 
of BD-22°3467 (Fig. The post-AGB evolutionary tracks 
were calculated by Miller Bertolami (priv. comm.) labeled with 
the stellar mass (in Mq). 

It is worthwhile to mention that there is no agreement with 
the HIP PARCOS distance a t the m aximum value of log § = 7.83 
given bv iHerald & Bianchil (l2002h either. The calculated spec- 
troscopic distance (using M - 0.68 Mq, R - 0.017 R©, and 
the respective model flux) is then D - 226 pc, about twice 
the HIPPARCOS distance. Although trigonometric parallaxes in 
general have to be regarded as hard constraints, orbital move- 
ment in binary systems such as BD-22°3467 may bear on 
the precision of HIPPARCOS measurements (cf. Rauch 20C)3, 
for LB 3459 = TYC 9166-00716-1). Figures |6] and [T] clearly 
demonstrate that the spectra cannot be reproduced with a grav- 
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Table 3. Mass, luminosity, radius, distance, and height above the 
Galactic plane of BD-22°3467. 



post-AGB evolution post-EHB evolution 



M/Mo 


0.52 ± 0.05 


0.48 ± 0.05 


log L/Lo 


1 50+"" 


1 00+" ^ 

^•""-0.80 




030+° '"^ 

"■"-'"-0.009 


029+" "'^ 

"•"^=^-0.009 


O/pc 




361+"^ 


///pc 


233!i 


224+^' 



ity of log g — 1.1 or higher, while a good fit is achieved for 
log ^ = 7.2. The upcoming GAIA mission will provide a highly 
precise parallax for BD-22°3467 which will be able to rule out 
the uncertainties. 



5. Results and discussion 

We calculated a grid of TMAP NLTE models atmospheres 
for BD-22°3467, the exciting star of the ionized nebula 
A 35. We determined Teff^SO+lOkK and logg^7.2±0.3 
and th us confirm Teft ^SOkK found by iHerald & Bianchil 
('2002'). Nevertheless, we found a somewhat lower log^ than 
[Herald & Bianchi ( 2002, log g-1.1). For the first time in this ob- 
ject, abundances from detailed line-profile fits were determined 
for Ar, Cr, Mn, Co, and Ni (Tab.©. For F Ne, Na, Mg, Al, Ca, 
Sc, Ti, and V, upper limits could be derived. 

The element abundance pattern (Fig. [T6] l is probably the 
result of the interplay of radiative levitation and gravitational 
settling. Compared to t hose objects in the DAO sample from 
lGoodetal.1 (12004 12005 ') that have the same T^fi and logg within 
error limits, BD-22°3467 has a higher He an d a lower Si abun- 
dance. Cr, Mn, Co, and Ni (not measured by iGood et ani2005h 
are strongly overabundant, likely due to the force of the radia- 
tion field. Radiative levitation may be a reason but the question 
remains why Fe is about solar Detailed diffusion calculations 
that consider these elements are highly desirable investigating 
this p henomenon. A com parison with existing diffusion calcula- 
tions (IChaver et al.lll995a.b . Fig.[T6]l shows (within error limits) 
that the metal abundances partially agree with predictions for 
DA-type WDs. Only N, Mg, S, and Ca are more than one dex 
lower than predicted. 

Until recently, BD-22°3467 was believed to be the central 
star of a planetary nebula and, thus, a post-AGB evolution of an 
intermediate-mass star (0.8 Mp ~ M j nitiai - 8.0 Mq) appeared to be 
a matter of course. iFrew & Parke r (2010) showed that the visi- 
ble, PN-like nebula may be a combination of shock excitation 
and photoionization of the ambient interstellar gas. 

BD-22°3467 provides enough high-energy pho- 
tons, (4.4 ■ 10^^ photons /s/cm^ at energies >13.6eV and 
6.4 ■ lO^^photons/s/cm^ at energies >54eV calculated from 
our final model which give totals of A?i3 = 2.3 ■ lO'*^ photons /s 
and A^54 = 3.3 ■ 10"*^ photons /s, respectively) to ionize 
the surrounding interstellar gas (cf. iRauch et"an l2004 l). For 
N54/N13 < 0.01, no Hen emission from the nebula can be ex- 
pected. If we assume an electron temperature of = lOOOOK 
in the nebula, negligible extinction, d - 361 pc (Table |4|i, 
and an average angular diameter of 750", we derive that 
BD-22°3467 has a Stromgren sphere with a density of 16 cm"^, 
log(FH/3/erg/s/cm^) = -11.23, and an ionized nebular mass of 
0.5 Mq. The absolute Hp fl ux agrees very well with the value 
given bv lAcker et aO (1 19921 log Fup ^ -U .3). 



Table 4. Results of this work compared with lHerald & Bianchil 
(I2002h . whose He abundance we adopted. 







this work 


H & B 




NLTE code 


TMAP 


TLUSTY 


primary 


7-eff/kK 


80± 10 


80 ±3 




log 8 


7.2 ±0.3 


7 7+0.13 

-0.18 




-'*-hle\-'*-0/ 


0.47 


0.3 - 0.5 




Xr(Xr,) 


1.2 X 10"' 


5.0 X 10"' - 10"' 




/V \T I /V /~i 1 

-'»-JN V-' *-0^ 


1.8 X 10"^ 


< 10"' 




^^OV^*^©.' 


2.5 X 10"' 


5.0 X 10"' - 10"^ 






~1 






Xx, fx 'I 


< \ 






Xx. fX ^ 


< \ 






x.« ex "1 


< \ 






X»,fX ^ 
^AlV,^©/ 


< \ 








7.5 X 10"' 


< 10"2 






^ 1 A V 1 0^3 








< 1 .Z X iU 






A. At (.A.©; 


1 9 






AcalAoj 


< i 






Asc(,Ao; 








ativa©; 


^ 1 






Av(,Aq; 


V. 1 






Xcr(Xo) 


70 






XMn(X©) 


35 






XFe(Xo) 


1 


0.1-0.5 




-^CoV-^©/ 


150 






XNi(Xo) 


5 






£^B-V 


0.02 ± 0.02 


0.04 ±0.01 




log«HI 


20.7 ± 0.1 


20.9 ± 0.1 






361+"^ 


163 (HIPPARCOS) 




M'' / Mo 


0.48 ± 0.05 


"■-'-0.4 




^"/R© 


0.029!»;°i^ 


0165+° '" 

"•"^"-'-0.006 




i?70"/Ro/kpc 


0.080!«;°°t 


0.1 ±0.01 




log L"/Lo 


1 00+"-^" 

^•""-0.80 


1.0 ±0.4 


secondary 


T-eff/kK 


5 


5 




log g 


3.5 


3.5 



our values for post-EHB evolution 



The neb ular spectrum by Acker & Stenholm obtained for 
the catalog (Acke r et"an 1 1 9921) is weak and, consequently, the 
Fua/Fuji ratio is overestimated, and the derived extinction of 
c - 0.7 is very uncertain. F^pjjjj^^Q^j ^/Fu0 - 2 indicates 

low excitation. Both N 11 /l/l 6548, 6584 A are visible and, with 
standard assumptions for the nebula but without any correc- 
tion for unseen higher ionization stages of nitrogen, its nebu- 
lar abundance may be estimated to be 12 + logN/H - 1.1. 
The Si II /i 6716 A / Si 11 /1 6731 A line intensity ratio is at the 
low-density limit and we can conclude that the nebular electron 
density is < 100 cm-'. In addition, we may estimate that the 
nebular oxygen and sulfur abundances are 7.9 and 6.4, respec- 
tively. It is worthwhile to note that the ionized mass of 0.5 Mq 
is of about the value usually assumed for the average FN mass 
of 0.2 - 0.3 Mq. The reason is simply the fairly low luminosity 
of BD-22°3467. Such a low-luminosity ionizing star makes a 
Stromgren sphere that may well mimic a genuine PN. It is thus 
possible that similar objects have been misclassified as PNe. 

Because the time since its departure from the AGB is rel- 
atively long for BD-22°3467 due to its low mass (Tab. O, 
all matter ejected previously on the AGB (as well as the PN) 
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Fig. 16. Top: Photosphe ric abundances of B D-22°3467 com- 
pared with solar valu es jAsplund et al. l 2009t) . and mean values 
from three objects of iGood et alj (120051 blue triangles) that are 
located within the error ellipse (cf. Fig. [15] green squares there). 
Arrows indicate upper limits. Bottom: Comparison of number 
ratios compared with predictions of diffusion calculations for 
hydrogen-rich (DA-, b lue squares) and he lium-rich (DO-type, 
green triangles) WDs dChaver et all 1 1995 alibi) with = 80kK 
and log § = 7.2. 



m ay have dissipa t ed int o the ambient ISM. Thus, the result 
of iFrew & Parkeii (1201 Ol) does not stringently exclude a post- 
AGB evolu tion. However, an a lternative stellar evolution may 
be possible. Berge ron et al.l (ll994) discussed that the progenitors 
of most DAO-type WDs are pos t-EHB stai's (M k 0.48 Mq). 
Recently, iGianninas et al.l (1201 Ol) presented a spectral analysis 
of 29 DAO- and 18 hot DA-type WDs and found that their 
DAOs (Fig. ITSl l are hotter and more massive that previously 
determined. Accordingly, they concluded that most of these 
( "a rather mixed bag of objects ") are products of post-AGB 
evolution. BD-22°3467 is located amongst the low-mass stars 
of their sample (Fig. ITTl l and consequently, a post-EHB evo- 
lution cannot be ruled out. In this scenario, the stars evolve 
directly from the EHB t o the WD state (AGB -manque stars, 
iGre^pio & Renzini|[T990h . iBergeron et al.1 ( 11994.) proposed that 
a weak r nass-loss may keep enough he lium in the line-forming 
regions (Unglaub & Bues ' 19981 l200d b to exhibit strong He ii 
absorption lines. Figure [T7] shows BD- 22°3467 in com pari- 
son to post-EHB evolutionary tracks (Dor man et alJll993l) . One 
can estimate that an M a; 0.48 Mq track (with solar helium 
content on the horizontal branch, HB) will match the posi- 
tion of BD-22°3467. (Dorman's tracks for masses higher than 
> 0.475 Mq and solar helium content are not fully calculated 
down to lower T^s and higher log g.) On its way from the EHB 
to its present position, i.e., towards higher T^^ and higher \ogg, 
gravitational settling, radiative levitation, and a weak wind are 
interacting and, thus, the atmospheric abundances may deviate 
from the HB abundances. 
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Fig. 17. Location of BD-22°3467 in the Teff-log^ diagram 
compared with sd(0)B-type s tars cl ose to the EHB taken from 
lEdelmanri| (blue, open squares l2003h and DAO-type WDs taken 
from Gianninas et al.l (red diamonds and green, filled squares 
'20 l ot cf. Fig. lTsTl. T he post-EHB evolutionary tracks are from 
.Dormanet al.lIl993L Y = 0.288 ^ Yq) and labeled with the stel- 
lar mass (in Mq). 
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Table 1. Log of the UV observations of BD-22°3467. The 
FUSE observation was performed with the LWRS (low- 
resolution aperture, resolving power/? = A/AA ^ 20000), STIS 
with grating E140M (R k 45 800), and lUE in low-resolution 
mode (R ~ 300). 



Instrument 


Obsid 


Start time 


Exp. time 






(UT) 


(s) 


FUSE 


P1330101000 


2000-05-20 20:27:37 


4416 


STIS 


O4GT02010 


1999-04-17 21:14:49 


2050 




O4GT02020 


1999-04-17 22:37:03 


2 800 




O4GT02030 


1999-04-18 00:16:10 


2740 




O4GT02040 


1999-04-18 01:52:54 


2 740 


lUE 


SWP44121LL 


1992-03-06 04:50:23 


1440 




LWP22518LL 


1992-03-06 05:26:54 


1440 



Table 2. Statistics of H - Ar model atoms used in our calcula- 
tions 



levels levels 



ion 


NLTE 


LTE 


lines 


ion 


NLTE 


LTE 


lines 


Hi 


12 


4 


66 


Navi 


43 


10 


130 


Hii 


1 







Navil 


1 








He: 


5 


98 


3 


Mg III 


1 


34 





He II 


16 


16 


120 


Mg IV 


31 





93 


Hem 


1 







Mg v 


15 


37 


18 


Cm 


6 


61 


12 


Mg VI 


1 








Civ 


54 


4 


295 


Aim 


1 


6 





Cv 


1 








Aliv 


15 


2 





Nlll 


1 


65 





Alv 


1 


16 





Niv 


16 


78 


30 


Alvi 


14 


24 


16 


Nv 


54 


8 


297 


Al VII 


1 








Nvi 


1 








Sim 


3 


31 


1 


Oiii 


3 


69 





Si IV 


16 


7 


44 


Oiv 


18 


67 


39 


Siv 


25 





59 


Ov 


90 


36 


610 


Si VI 


1 








Ovi 


54 


8 


291 


Fill 


1 


9 





Ovii 


1 








Piv 


15 


36 


9 


Fill 


1 


6 





Pv 


18 


7 


12 


Fiv 


1 


10 





Pvi 


1 








Fv 


15 


91 


31 


Siv 


6 


94 


4 


Fvi 


12 


115 


16 


Sv 


21 


89 


37 


Fvii 


1 








Svi 


18 


19 


48 


Neil 


1 


33 





Svii 


1 








Neiii 


3 


43 





Ariv 


1 


349 





Neiv 


3 


37 





Ar V 


32 


329 


38 


Ne V 


20 


74 


35 


Arvi 


16 


168 


21 


Ne VI 


1 








Arvil 


40 


112 


130 


Nam 


1 


186 





Ar VIII 


1 








Naiv 


1 


237 













Nav 


8 


42 


9 


total 


741 


2837 


2514 



Table 2. Continued for Ca - Ni. The so-called "sample lines" 
are combined to "super lines" by IrOnIc. 



ion 


super levels 


super lines 


sample lines 


Caiv 


6 


16 


20291 


Cav 


6 


21 


141956 


Cavi 


6 


19 


114 545 


Cavil 


6 


21 


7 1608 


Ca VIII 


6 


20 


9 124 


Calx 


1 








Sc IV 


6 


20 


15 024 


Scv 


6 


21 


261 235 


Scvi 


6 


19 


237 271 


Scvii 


6 


20 


176 143 


Sc VIII 


6 


21 


91935 


Scix 


1 








Tiiv 


6 


19 


1 000 


Tiv 


6 


20 


26654 


Tivi 


6 


19 


95 448 


Tivil 


6 


20 


230 618 


Tivm 


6 


21 


182 699 


Tiix 


1 








Viv 


6 


19 


37130 


Vv 


6 


20 


2 123 


Vvi 


6 


19 


35 251 


Vvii 


6 


19 


112 883 


V VIII 


6 


20 


345 089 


Vix 


1 








Criv 


6 


20 


234 170 


Crv 


6 


20 


43 860 


Crvi 


6 


20 


4406 


Cr VII 


6 


19 


37 070 


Cr vm 


6 


20 


132 221 


Crix 


1 








Mniv 


6 


20 


719387 


Mn V 


6 


20 


285 376 


Mn VI 


6 


20 


70116 


Mn VII 


6 


20 


8 277 


Mn VIII 


6 


20 


37168 


Mnix 


1 








Feiv 


6 


20 


3 102 371 


Fev 


6 


20 


3 266 247 


Fe VI 


6 


20 


991 935 


Fe VII 


6 


20 


200455 


Fe vm 


6 


18 


19 587 


Feix 


1 








Co IV 


6 


20 


552 916 


Cov 


6 


20 


1469 717 


Co VI 


6 


18 


898 484 


Co VII 


6 


19 


492 913 


Co VIII 


6 


20 


88 548 


Co IX 


1 








Niiv 


6 


20 


2512561 


Niv 


6 


20 


2766 664 


Ni VI 


6 


18 


7 408 657 


Ni VII 


6 


18 


4195 381 


Ni VIII 


6 


20 


1 473 122 


Niix 


1 








total 


279 


884 


33 219 636 




Fig. 2. Temperature and electron density stratification along with ionization fractions of all elements in our final model. 



